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Abstract 



The cosmological constant problem is studied in a two component cosmological 
CJ i model. The universe contains a cosmological constant of an arbitrary size and 

sign and an additional component with an inhomogeneous equation of state. It is 
shown that, in a proper parameter regime, the expansion of the universe with a 
large absolute value of the cosmological constant may asymptotically tend to de 
Sitter space corresponding to a small effective positive cosmological constant. It 
^ ■ is argued that such a behavior can be regarded as a solution of the cosmological 

^N , constant problem in this model. The mechanism behind the relaxation of the 

vQ I cosmological constant is discussed. A connection with modified gravity theories 

cn ' is discussed and an example of a possible realization of the cosmological constant 

t^^ ■ relaxation in f{R) modified gravity is described. 
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1 Introduction 

The state of accelerated expansion of the present Universe seems to be better and 
better confirmed by the cosmological observations [H El [3] . The question of the 
dynamical mechanism responsible for the accelerated expansion, however, still lacks 
its definite answer. The last decade has witnessed the arrival of numerous models 
of the accelerated expansion including dark energy, braneworld models, modified 
gravity and many others [3]. An interesting fact is that when confronted against 
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the observational data, a simple ACDM model, in which the cosmological constant 
(CC) is the cause of acceleration, fits the data very well. When the conceptual 
simplicity of the ACDM model is taken into account, it is easy to understand why 
it is a benchmark model for the analysis of cosmological observations. While it is 
quite clear why the ACDM model is so appealing from the observational side, its 
status from the fundamental theoretical perspective is much more problematic. It 
is by now a notorious fact that the CC value predicted in Quantum Field Theory 
(QFT) differs from the observed value by an embarrassingly large number of orders 
of magnitude [5] . The problem of explaining the observed value of the cosmological 
constant is therefore one of the largest challenges in theoretical physics O El [7] . 
This problem is sometimes refered to as the old cosmological constant problem. It 
is further exacerbated by the fact that in all other approaches to the problem of 
the accelerated expansion of the universe it is assumed that the CC problem is 
somehow solved. 

Many attempts to solve the CC problem have been made during several last 
decades El E]- However, so far none of them has provided a fully satisfactory 
solution of the CC enigma. The most frequent problem that models of various 
sorts encounter is the necessity of fine-tuning. The parameters of the model have 
to be chosen with extraordinary precision in order for the model to lead to the 
resolution of the CC problem. Even a very small deviation from these fine-tuned 
values disrupts the efficiency of the proposed mechanisms. 

In this paper we propose a dynamical cosmological model with a specific regime 
in which it is possible to contemplate the resolution of the CC problem. The 
proposed model is simplified insofar that it does not contain all the (matter or 
radiation) components that naturally participate in the evolution of the universe. 
However, given the difficulty of the CC problem and its resilience to different 
attempts of solution, it seems preferable to first concentrate on the very mechanism 
which could produce the observed value of the effective cosmological constant for 
a universe with values of A comparable to those predicted in QFT. 

There exists additional problem related to the size of the cosmological constant 
energy density (or more generally the present dark energy density). Namely, the 
observational data reveal that the energy density of matter (which at present epoch 
comes predominantly from nonrelativistic matter) is of the same order of magnitude 
as the CC (present DE) energy density. Dark energy and nonrelativistic matter 
scale differently with the expansion of the universe and it is quite remarkable that 
at present epoch these two energy densities are comparable. This problem is also 
called the cosmic coincidence problem. In this paper we are primarily concerned 
with the problem of the size of the cosmological constant, whereas the cosmic 
coincidence problem is not addressed. 

The principal aim of this paper is to study a two component cosmological model 
which possesses a cosmological constant of a large absolute value. We investigate 
conditions under which the asymptotic expansion in this model is of de Sitter type 
where the asymptotic value of the Hubble function corresponds to a small value 
of the effective cosmological constant {H'^sym — ^e///3). In the studied models 
the solution of the cosmological constant is understood as a situation in which the 
universe ends up in an asymptotic de Sitter regime at large scale factor values 
characterized by a small cosmological constant (small in a sense that Ag// <^ |A|). 



The absence of fine-tuning in the model studied in this paper is achieved if the 
parameters of the model do not have to cancel to many decimal places to lead to 
the solution of the CC problem. The dynamical process in which the universe in 
the studied model tends to de Sitter space with a small Ae// is also referred to as 
the relaxation of the cosmological constant. In the following sections we present 
a two component model in which the relaxation of the cosmological constant is 
realized. 



2 The model set-up 



We consider the cosmological model described by the FRW metrics containing 
two components: the cosmological constant with the energy density /ja and an 
additional cosmological component with the energy density p. Throughout this 
paper the universe is assumed to be spatially flat (/c = 0). The Friedmann equation 
for this cosmological model is 

H^ = ^{PK + P). (1) 

The equation of state (EOS) of the cosmological constant is standard, pA = —p\. 
On the other hand, for the equation of state of the second component we take 

p = wp- 3Coii'"+' , (2) 

where p is the pressure of the second component and w, Co and a are real parame- 
ters. In the considerations given in this paper we limit ourselves to positive values 
of the parameter Co- Both components satisfy the equation of continuity which for 
the cosmological constant results in pA = const and for the second components 
reads as 

d/> = -3(p + p)-, (3) 

a 

where a denotes the scale factor of the universe. 

Although the preceding equations of this section describe a simple cosmologi- 
cal model with interesting dynamical regimes, as it will be shown in the following 
section, it is clear that in Eq. ^ lies the nonstandard content of this model which 
requires physical motivation. This equation of state can be described as an in- 
homogeneous equation of state in the framework of Ref. [8] (the same authors 
consider very similar EOS in a different context). An interesting example where 
an inhomogeneous dark energy EOS is relevant in the structure formation process 
is given in [9]. Our primary aim is to demonstrate a mechanism of this model 
which allows a universe with a large |A| to end up in a de Sitter regime with a 
small positive Ag//. An important issue of a more fundamental basis leading to ([2]) 
is left for future work. Still, an illustrative example of a possible realization of the 
relaxation mechanism in terms of f{R) gravity is given in section [H In the remain- 
der of this section we outline two physical frameworks [8] which give motivation for 
the inhomogeneous equation of state of the type ^. It is important to stress that 
our approach in the study of the model is mainly phenomenological. We primarily 
focus on the mechanism of the CC relaxation whereas the model studied in this 



paper is mainly considered as a framework in which the said machanism could be 
studied. 

The form of Eq. ^ is intentionally chosen to emphasize similarity with bulk 
viscosity. Indeed, the identification C = Co-f^" brings ([2]) in the form p = wp — 2>HC, 
which is the standard form for the bulk viscosity eff'ects of imperfect cosmological 
fluid in FRW universe [lOtlllj. However, it should be noted that the dependence of 
C, on H, which is not a state variable of the fluid, does not correspond to standard 
bulk viscosity. For the value a = we recover the bulk viscous imperfect fluid 
with the constant coefficient ("• Therefore, we are motivated by the bulk viscosity, 
we consider its generalization and then proceed phenomenologically. Here it would 
be preferable to call the inhomogeneous term the nonlinear (bulk) viscosity |12j . 
Therefore, a possible identification of the second component might be as an im- 
perfect cosmological fiuid with nonlinear viscosity having a power-law dependence 
on the Hubble parameter H. 

An alternative view is to interpret ([2]) as an inhomogeneous equation of state 
coming from modified gravity or braneworld models (see Appendix in [H]). The 
considerations given in section H] lend support to this interpretation. 

In the following section we focus on the dynamical regimes of the system of 
equations ([l|), ([2]) and ([3]) and the related phenomena including the possibility of 
the relaxation of a large cosmological constant (large in terms if its absolute value) . 

3 The model dynamics 

The model defined in the preceding section is now analyzed in detail. Combining 
([l]) and ([2]) with ([3]) we obtain a dynamical equation for the Hubble parameter 

dH^ + 3(1 + ^)^ (h^ - 5^^^ - ^!^(i/2)(a-M)/2\ ^^ (4) 

a \ 3 \ + w J 

The analysis of this equation is further simplified by the introduction of the fol- 
lowing notation: 

h = {H/Hxf, s = a/ax, X = SttGpa/SH],, ^ = SttGCqH'^-^ /{I + w) . (5) 

Here Hx denotes the value of the Hubble parameter at, in principle arbitrary, 
value of the scale factor ax- Let us further stress that although we assume Co > 0, 
the parameter ^ may take values of both signs if we also allow the values u^ < — 1. 
Applying the described change of notation, Eq. ^ acquires the form 

s^ + 3{l + w){h-X-^h'^) = 0, (6) 

as 

with the initial condition h{l) = 1. 

The inspection of Eq. ([6]) reveals that the value of the parameter a may sig- 
nificantly infiuence the type of dynamics of the Hubble parameter. The values 
a = —1 and a = 1 are specific points at which we expect the change of dy- 
namical behavior. Therefore we analyze five characteristic intervals/points for a: 
(— oo, —1), —1, (—1, 1), 1 and (l,co). For each of the intervals/points we make an 
analytical treatment at one value of the parameter a and, when necessary, support 
it with numerical calculations. 



3.1 a < —1: the relaxation mechanism for a large cos- 
mological constant 

The analysis of the interval a < —1 reveals a dynamical mechanism for the re- 
laxation of a large cosmological constant to a much smaller effective CC value. 
The contents of this subsection comprise the main results of the present paper. 
We start the analysis of the Hubble dynamics with an analytical treatment for a 
representative value a = —3. Eq. ([6]) now becomes 

hdh , ds 

-3{l + w)-. (7) 



The integration of the left hand side of Eq. ([7]) is determined by the zeros of its 
denominator which we denote by /i*i,2, i-e. h? — \h — £^ = {h — /i*i)(/i — /i*2)- Their 
specific values are 

h.i = \[\ + \l\'' + ^dj , (8) 

1 



/i*2=2(A-^A2 + 4ej . (9) 

Generally we have /i*i — /i*2 > 0. The solution of ([7]) is of the form 

h-K_,Y^(h-K2Y'^s-''~'^-\ (10) 



1 — /l*l / \ 1 — /i*2 

where Ai = /i*i/(/i*i — /i*2) and A2 = — /i*2/(^*i — h^,2)■ Let us further study 
separately cases of positive and negative values of the parameter ^. 

For w > —1 the value of the parameter ^ is positive. In this case we have 
/i*i > and /i=k2 < which leads to ^41 > and A2 > 0. Eq. (fTOl) readily provides 
information on the asymptotic behavior of the parameter h. Namely, for small 
values of the scale factor, s — > the function h diverges 

/i ~ (1 - /l,i)^i (1 - /l,2)^^S~'^'+"^ . (11) 

For large values of the scale factor the function h tends to a constant value 

lim h = h^:i . (12) 

As stated in (I12p . the parameter h (equivalently the Hubble parameter squared 
H'^) asimptotically tends to a constant value at large values of the scale factor. Let 
us consider the case when A is negative and very large in absolute value compared 
to ^/^, or, more precisely, A^ » 4^. The square root in the expression ([8]) can now 
be expanded and we obtain 

or, equivalently, 



^*i - Mj ' (13) 



2 ^ 24vrGCo ^ 3Co 

*i {l + w)\A\ ' {1 + w)\pa\- ^ ^ 

For p\ very large in absolute value and negative, the universe asymptotically 
tends to a small value of H^ which can be interpreted in a straightforward manner 
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Figure 1: The evolution of h = H^ / H\ as a function of the scale factor for different 
values of the exponent a. The value of a strongly influences the asymptotic value of 
h at large a, whereas the behavior at small a is not affected by a. The values of the 
parameters used are A = —1000, ^ = 0.01 and w = —0.9. 



as a small value of the effective cosmological constant. It is sufficient that A is neg- 
ative and large in absolute value and that £, is sufficiently small. The model under 
study in this paper, therefore, provides a dynamical mechanism for the relaxation 
of a negative A with a very large absolute value. The dynamical mechanism does 
not incorporate fine-tuning of model parameters and it can be considered as a 
solution of the cosmological constant problem in this model for a negative large A. 

Although the choice a = —3 gives an analytically tractable example of the re- 
laxation of a negative A, we further support the findings of the preceding paragraph 
with numerical solutions of Eq. ^ for other values in the interval a < — 1. 

In Fig. [T]we present the evolution of the variable h = H"^ /H'j^ as a function of 
the scale factor for different values of the exponent a. The most striking feature of 
the dynamics of h for all studied values of a, is the abrupt transition between two 
asymptotic regimes. This feature could not be properly addressed from the study 
of the asymptotic regimes alone, but it requires a numerical treatment to be fully 
appreciated. The value of exponent a does not affect the asymptotic evolution at 
small a, but it crucially affects the large a asymptotic behavior. 

Figure [2] depicts the dependence of the dynamics of h on the CC parameter 
A. It is evident that the value of A affects both asymptotic regimes (at large and 
small scale factor values) as well as the onset of the abrupt transition between two 
regimes. 

In Fig. [3] we study the dependence of the behavior of h on the parameter ^. 
From the figure it is clear that the dynamics of h at small values of scale factor 
does not depend on ^, but at large values of a the asymptotic value of h is strongly 
influenced by the value of £,. 

Finally, in Fig. H] we present the dynamics of /i as a function of the scale factor 
a for various values of the parameter w. The plots in the figure reveal that the 
behavior at small a is strongly affected by w, whereas the asymptotic behavior at 
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Figure 2: The dynamics of h 
values of the CC parameter A. The value of A affects the behavior at small values of a, 
large values of a and the onset of the transition between two regimes. The values of the 
parameters used are a = —3, S, = 0.01 and w = —0.9. 
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Figure 3: The dependence of /i as a function of the scale factor on the parameter ^. It is 
evident that the behavior at small values of a does not depend on ^, but the asymptotic 
dynamics at large a is strongly affected by the value of the parameter ^. The values of 
the parameters used are A = —1000, a = —3 and w = —0.9. 
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Figure 4: The evolution of /i as a function of scale factor for different values of the 
parameter w. The dynamics at small a is clearly affected by the value of w and the 
dynamics at large values of a is insensitive to parameter w. The values of the parameters 
used are A = —1000, a = —3 and ^ = 0.01. 



large scale factor values does not depend on w. 

Next we return to our analytically tractable case of a = —3, but this time we 
consider a positive value of the cosmological constant A. We also choose w < —1 
so that the parameter ^ becomes negative. In this setting we have /i*i > and 
/i*2 > which results in ^i > and A2 < 0. The asymptotic behavior of the 
Hubble parameter now obeys the following laws: 



lim h 



h*i 



and 



lim h = hif2 



(15) 



(16) 



We see that for large values of the scale factor the scaled Hubble parameter tends 
to a constant value /i*2- For a very large values of A (such that A^ ^> —4^) the 
universe asymptotically acquires a small value 



or, equivalently. 



H: 



h*2 



247rGCo 



A' 



3Co 



*2 



(17) 



(18) 



{l+w)A {1 + w)pa' 

For a very large value p\ and a small value of |^| the asymptotic value (jlSp is very 
small. This small value of H^ can be directly interpreted as a small effective positive 
cosmological constant. As for the case of negative A, we have at hand a dynamical 
mechanism of the cosmological constant relaxation which does not incorporate fine- 
tuning of model parameters. Therefore, within the model studied in this paper, we 
have presented a solution of the CC problem for a positive cosmological constant. 
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Figure 5: The dependence of h on scale factor for different values of the exponent a. 
The graphs reveal that the behavior of h at small values of a is insensitive to the value 
of a, whereas at large values of a there is a strong dependence of the asymptotic value 
of h on a. The values of the other model parameters used are A = 1000, ^ = —0.01 and 
w; = -l.l. 



It is important to notice that p in this case must be negative. This fact implies 
that the component with an energy density p is an effective description of some 
other dynamical mechanism, possibly modified gravity. 

As for the case of negative CC, we further support the analytical treatment for 
a = —3 with numerical analyses for other values in the interval a < — 1 and other 
model parameters. 

In Fig. [5] the dependence of h on the scale factor for different values of the 
exponent a is depicted. The behavior of h for small values of the scale factor 
shows no dependence on the exponent a. At large values of the scale factor the 
asymptotic value of h is strongly influenced by the value of a. As in the case 
with the negative CC with a large absolute value, here the transition between 
dynamical regimes at small and large values of a is abrupt. Another significant 
difference compared to the case of negative CC is that both at small and large 
values of the scale factor the expansion is of de Sitter type (for negative CC we 
have a de Sitter-like expansion only at large values of the scale factor) . 

The dependence of h on the scale factor for various values of the CC parameter 
A is given in Fig. [6l The value of A affects the dynamics at small and large values 
of a and the onset of the abrupt transition between two regimes. 

The sensitivity of the dynamics of /i as a function of a on the parameter S, is 
presented in Fig. [71 The dynamics at small values of a is not affected by the 
value of ^, but at large values of the scale factor the asymptotic value of h depends 
strongly on ^. 

Finally, in Fig. [8] the dependence of h on the scale factor for several values of 
w is depicted. The plots in the figure reveal that the value of w does not affect 
the asymptotic values of h at small and large values of the scale factor, but they 
do affect the transition between these asymptotic values. For all studied values of 
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Figure 6: The dynamics of /i as a function of the scale factor for various values of the 
parameter A. The CC parameter A affects the behavior at small and large values of the 
scale factor as well as the onset of the abrupt transition between two asymptotic regimes. 
The values of the other parameters used are a = —3, ^ = —0.01 and w = —1.1. 
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Figure 7: The dependence of h on the scale factor for various values of the parameter ^. 
The dynamics at small a is not influenced by ^, but at large a the asymptotic behavior is 
sensitive to the value of ^. The values of the parameters used in this figure are a = — 3, 
A = 1000 and w = -1.1. 
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Figure 8: The dependence of /i as a function of a for different values of the parameter 
w. This parameter affects the form of transition between asymptotic values of h which 
are insensitive to w. The values of the parameters used are a = —3, A = 1000 and 
^ = -0.01. 

w this transition remains quite abrupt. 



3.2 a = -1 

The solution of dH) for a 



h 



-1 can be expressed in a familiar form 

A + e + (1 - A - 0^"^^^+"'^ • 



(19) 



This expression describes the universe with the cosmological constant energy den- 



sity ~ 8-^(A-|-^) and the matter component with the scaling ~ a' 



-3{l+w) 



Clearly, 



for w = our model fully mimics the ACDM model. If we wish to have a small 
effective value of the cosmological constant ~ A + ^, for a large value of A we need 
to have a cancellation of A and ^ which necessarily introduces fine-tuning. There- 
fore, for a = —1 we do not have an efficient mechanism for the relaxation of the 
cosmological constant. 



3.3 -1 < a < 1 

In this interval for the value of a we consider the case a = 0. The Eq. 

reads 

ydy 3 ds 

y^ -iy-X 2 s 

where y = h}'"^ . The solution of this equation can be cast in the form 



1 -y*] 



1 -y*2 



where 



y*i = ^{^ + \/e + 4A 



now 
(20) 

(21) 
(22) 
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1 Ic 


- ^/e + 4A j 


5 


(23) 


with Bi = y*i/(y*i - 


- y*2) and B2 = 


-y*2/{y*i - 


- y*2)- 


We again study two 


interesting cases 











1. We consider t;; < — 1 where we have ^ < 0. For the value of the cosmological 
constant A let us take a positive value which is sufficiently small, so that 
^^ » 4A. Now we have y=Ki > 0, y*2 < with y*i — y*2 > and B12 > 
0. The asymptotic behavior of the system is given by lims__»o y = y*i and 
lim^^ooy = 00 with y ~ (1 — y^.i)^^{l — 2/^,2)^^ s"^^^^'"-'''^- At large values of 
the scale factor we do not have a de Sitter regime although there is one at 
small values of a. 

2. Another interesting regime is obtained for w > —1 which implies ^ > 0. We 
study the case of a negative A which is again sufficiently small in absolute 
value i.e. ^"^ » 4|A|. Then we have y*i > 0, y*2 > with y^.i — y^,2 > and 
Bi > 0, B2 < 0. The asymptotic behavior of the system is lim^^oy = y*2 
and lim^-too y = y*i- From ([22]) it is easy to see that y*2 — — A/^ and y*i ~ ^. 
For a large value |A|, ^^ has to be even larger and therefore at large values of 
the scale factor we have H^ '^ y'^2 ~ ^'^ what does not correspond to a small 
value of effective CC. 

3.4 a = l 

For a = 1 and C 7^ 1 the solution of Eq. ([6]) acquires the form 

/. = ^ + (1 - ^) .-3(i+-)(^-«) . (24) 

There are two cases of interest for the CC problem: 

1. For w < —1 we have ^ < 0. Let us further consider the case A > with 
A/(l — ^) <C 1. Under these conditions the asymptotic behavior of h is the 
following: at small values of the scale factor we have liuis—^oh = A/(l — ^), 
whereas for large values of the scale factor we have lims^oo h = 00 with 
/i ~ (1 — A/(l — ^))s~3(i+u')(i-5)^ jj^ |.]^jg pg^gg ii^Qj-Q jg j-^Q asymptotic de Sitter 

solution for large values of a. 

2. Uw > -1 then ^ > 0. Let us take 1 - ^ < with ^ «; 1 and A/(l - C) «; 1. 
Then we again have lims^o ^ = A/(l — ^) and linis^oo ^ = co with h r^ 
(1 — A/(l — ^))s~3(i+u>){i~^)^ Again there is no realization of the scenario of 
interest. 

For a = 1 and ^ = 1 the solution for the dynamics of h becomes 

h = l + 3X{l + w)lns. (25) 

This solution does not asymptotically lead to de Sitter space at large scale factor 
values. 
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3.5 a > 1 

As a representative and analytically tractable case in this interval of the exponent 

a we study the dynamics of h for the value a = 3. The dynamical equation ([6]) 

can now be written 

dh , . ^ds , . 

Introducing the notation /i^ — h/S, + A/^ = {h — /i*i)(/i — /i*2) we have 

h*i = \{\ + \lii-'l^\. (27) 



which leads to the following solution: 

h — /i*i 1 — /i*i 
h — hif2 1 ~ /i*2 




^*i = i 7-t/^-^l ' (28) 



g3(l+«))C(/i*i-/i.2) _ og) 



Let us now restrict ourselves to the parameters satisfying 1 — 4^A > as a prereq- 
uisite for having real asymptotic values for h. Again we consider two cases: 

1. For w > —1 which implies ^ > the condition of reality of asymptotic values 
of h leads to the requirement A < 1/(4^). We have /i^,! — /i*2 > 0, /i*i > 
whereas /i*2 > for < A < 1/(4^) and K2 < for A < 0. For A > 
we have linis^o ^ = ^*i and lim^^oo ^ = /i*2- For A <^ 1/(^0 we obtain 
/i^,i ~ 1/^ and /i=k2 ~ A. This scenario is not fully satisfactory since A itself is 
the effective CC at large scale factor values and, therefore, A cannot have its 
natural QFT value. For A < there is no de Sitter behavior at large values 
of the scale factor. 

2. In the case when w < —1 we have ^ < 0. The condition of reality of 
asymptotic values results in a requirement A > 1/(4^). We further have 
/i*i — /i*2 > 0, /i*2 < and /i*i > for A > whereas /i*i < for 
1/(4^) < A < 0. The only case at which we could have a de Sitter behavior at 
large values of the scale factor is for A > 0. However, the asymptotic behavior 
is lims_>o h = /i*i and at large a there is no de Sitter behavior. 

4 The general conditions for a small effective 
CC and a hint from modified gravity 

Many results of this paper, and especially the main results presented in subsection 
13.11 stem from the analysis of the asymptotic behavior of the dynamics of h. The 
solutions for the function h{s) which tend to a small positive constant value at 
large values of the scale factor have been interpreted as solutions of the cosmological 
constant in our approach. In this paper we use a specific form of the inhomogeneous 
EOS for the component p with C{H) = CqH". A natural and important following 
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step would be to consider a broader class of functional behavior for C(^)- In general 
the dynamics of the Hubble function would then be governed by the equation 

rlh 

s^ +3{l + w){h-\- C'h^^^Cih)) = , (30) 

as 

where S,' is a constant. A prerequisite for a solution of the CC problem for the 
function C(^) is the existence of a small and positive root /i* <C |A| of the equation 

h-X-^'h^/^C{h) = 0. (31) 

An additional condition is that the function h asymptotes to /i* at large values of 
the scale factor. 

As already stated in section [2] and elaborated in [8], a theory behind the in- 
homogeneous EOS of the type ([2]) could be some formulation of modified gravity. 
In the remainder of this section we perform an analysis of a possible asymptotic 
behavior in a model of f{R) modified gravity and discuss its implications for the 
solution of the CC problem. For a review of the f{R) modified gravity see [13] and 



We consider a f{R) theory with a arbitrarily large cosmological constant energy 
density. In a universe with the FRW metric the dynamics of H is given by the 
equation 



3f{R)H^ - ^(Rf'iR) - f{R)) + 3HRf"{R) = SvtGm ^ A , (32) 

where prime denotes the differentiation of f{R) with respect to its argument and 
R = 12ij2 + 6^. Next we choose 

/(i?) = i? + 6i?2_^i__, (33) 

which satisfies the requirements of stability and positivity of the effective gravi- 
tational coupling |14t I15j . Here 6 > 0, /x and n are the parameters of the model. 
In general in this model, as well as in many other models of f{R) gravity, there 
is not flat solution (with R = 0). Furthermore, modified gravity theories can be 
subjected to stringent local gravity tests, e.g. measurements in the Solar system. 
We assume that the values of parameters n, b and /i used here are consistent with 
the bounds from local gravity tests. Next we focus on the asymptotic behavior of 
H and search for constant H solutions of Eq. (|32|) . We neglect all time derivatives 
and (132]) can now be written as 

Finally, for illustration purposes, we choose n = 1 and a possible asymptotic value 
of H is determined by the equation 

3 48 ' ^ ^ 
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the solutions of which are 

Hk2 = -J'i±\l{i] +TTT1- (36) 




For a negative A with a large absolute value (so that 3;U^/4A^ *^ 1) the + solution 
becomes 

The comparison with the results of subsection 13.11 shows a striking similarity with 
our model containing a component with inhomogeneous EOS. In the modified 
gravity model (I33p there is an asymptotic de Sitter behavior corresponding to a 
small Ae//. This finding strongly supports a conjecture that the relaxation of 
a large |A| is also feasible directly in modified gravity theories. The details of 
the relaxation mechanism in f{R) modified gravity theories will be elaborated 
elsewhere 1161. 



5 Discussion 

The preceding sections define our model, explain its inherent mechanism for the 
solution of the cosmological constant problem and outline its potential connection 
with the modified gravity theories. In this section we further discuss the aspects 
of the model which are of relevance for cosmological issues. 

There are at least two sources of motivation for the definition of component p 
in terms of an inhomogeneous EOS. The first one is nonlinear (bulk) viscosity and 
the second one is modified gravity. Although both of these possibilities contain 
many specific variants, in this paper the focus has been on the approach valid for 
both of them. 

As described in subsection 13.11 in suitable parameter regimes it is possible to 
have the solution of the CC problem for both positive and negative CC of large ab- 
solute value. Still, for negative A the energy density p remains positive throughout 
the evolution of the universe, whereas for positive A the energy density p must be 
negative. This is a strong signal that, at least for positive A, we should consider the 
component p as an effective description of a more fundamental dynamics, possibly 
coming from the modifications of gravity. 

In our treatment we assumed that Co > in order to stay aligned with a 
possible interpretation that inhomogeneous term in ([2]) might come from some 
sort of nonlinear (bulk) viscosity. This assumption immediately relates the signs 
of parameters 1 + w and ^. With this approximation, the relaxation of a large 
positive CC requires w < —1, i.e. the parameter of the EOS should be of the 
phantom type [l^. Since in modified gravity models it is possible to obtain the 
phantom- like effective dark energy [TH] , the interpretation of the component p as 
an effective description of the modification of gravity gains further support. The 
requirement Co > allows us to simultaneously treat both sources of motivation for 
the inhomogeneous EOS of the component p. The relaxation of this requirement, 
i.e. allowing for negative values of Co opens up space for scenarios that might be 
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realized in specific theories of modified gravity. The study of models with negative 
(^0 is an important challenge of future work. 

Quantum field theory provides various positive and negative contributions to 
A. Negative contributions come from the zero point energy of the fermionic degrees 
of freedom and some condensates. On the other hand, positive contributions come 
from the zero point energy of the bosonic degrees of freedom. In the analysis of 
the cosmological constant problem the size and interplay of positive and negative 
contributions is of considerable importance and so is the sign of the resulting A. An 
important conclusion of the present paper is that for any sign of A in the model with 
the inhomogeneous EOS ([2]) we can end up in a universe with a small and positive 
effective cosmological constant. Whereas our model provides mechanism for the 
relaxation of both positive and negative CC, the interpretation of the component 
p might differ. The present analysis shows that for the relaxation of a negative A 
with a large absolute value p could equally play a role of a real cosmic fluid with 
nonlinear viscosity or be an effective description of the modified gravity effects. In 
the case of large positive A the effective nature of p seems more plausible. 

The analysis of subsection 13.11 further reveals that for both signs of the cosmo- 
logical constant we can have two phases of accelerated expansion connected with an 
abrupt transition between them|^. It is important to stress that these two phases 
of the accelerated expansion correspond to very different energy densities. It is an 
open question for future research whether these two phases of accelerated expan- 
sion and the abrupt transition between them could be useful in the description of 
the inflationary dynamics and the graceful exit. 

From the analysis of the subsection 13.11 we can see that the asymptotic value of 
the Hubble function H depends on parameters ^ and A i.e. on their ratio. Let us 
discuss our expectations of the size of these parameters. In the remainder of this 
paragraph we use the terms "large" and "small" loosely for illustration purposes. 
For the parameter ^ we do not expect to be "large". Namely, it describes either 
effects of nonlinear viscosity or deviations from general relativity. A natural size 
of £, could be described as "small" . On the other hand, from QFT we expect the 
absolute value of A to be "large" . When ^ is "small" and A is "large" , the resulting 
value of H^ is very small especially compared to | A|/3. Therefore, for proper values 
of a and for the parameters ^ and A taking their "natural" values we have a small 
asymptotic value of H^. Within our model this solves the cosmological constant 
problem for any sign of A. 

One of the most important issues towards a complete cosmological model which 
would incorporate the mechanism of the CC relaxation is the addition of matter 
and radiation components. The cosmological model should reproduce the eras of 
radiation domination and matter domination to be consistent with the available 
observational data. A dedicated analysis is required to account for details of the 
model with matter and radiation components. It will be particularly interesting 
to see how the addition of matter and radiation components interacts with the 
abrupt transition characteristic for the CC relaxation mechanism. The cosmic 
coincidence problem could be possibly addressed only in such a full cosmological 



^For a negative A we need to have w < — 1/3 to have an accelerated expansion at small values of the 
scale factor. 
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model. Although the full treatment is needed for the understanding of the entire 
dynamics of the full model, it is possible to argue that the addition of matter 
and radiation components will not affect the asymptotic behavior of the model 
and therefore the very CC relaxation mechanism. Namely, the energy densities 
of matter and radiation components decrease quickly with the expansion and at a 
sufficiently large scale factor value they become negligible. Then we are effectively 
back to the two component model presented in this paper and asymptotically we 
have the CC relaxation. 

Finally, the focus of this paper is on the verification of the very effect of the 
relaxation of the large CC. There are many important questions that need to be 
addressed before our model with the inhomogeneous EOS could become a complete 
cosmological model. As already stated, it is important to learn if and how the 
behavior of the model changes when other cosmological components like radiation 
or nonrelativistic matter are added. Another question of considerable importance 
is which types of the inhomogeneous EOS apart from the one studied in this paper 
are capable of the relaxation of the large CC. The connection with the modified 
gravity theories seems especially worth pursuing. The growth and stability to 
perturbations is a relevant question too and so is the realization of the relaxation 
mechanism in astrophysical gravitationally bound systems. These questions are 
left for future work. 

6 Summary and conclusions 

The cosmological constant problem is a spot in theoretical physics landscape where 
the inadequacy of standard theoretical approaches is evident. An unconventional 
new ingredient is clearly called for. It is unclear, however, how big a deviation 
from the standard principles this new ingredient should represent. In this paper 
we have presented a simple approach based on a cosmological component with an 
inhomogeneous equation of state. The new ingredient is an inhomogeneous term 
in ([2]) which can be interpreted as nonlinear viscosity or the effect of modified 
gravity. In a particular parameter regime, which by itself requires no specific fine- 
tuning, the universe with a cosmological constant of any sign and an arbitrarily 
large absolute value asymptotically ends up in de Sitter regime with a small value 
of the effective cosmological constant. This result provides a solution of the CC 
problem without the need of fine-tuning. A preliminary analysis of a f{R) modified 
gravity theory lends support to the claim that the mechanism of the CC relaxation 
studied for a model with EOS ^ also functions for modified gravity theories. 
The main results of this paper exemplify a scenario in which a large A coming 
naturally from QFT calculations coexists with a small asymptotic value of the 
effective cosmological constant. The said results further open a possibility that the 
measured value of the cosmological constant is not the value coming from QFT, 
but it is determined by the QFT value. Apart from the intrinsic value of these 
results, they also allow us a bit different perspective on cosmological parameter 
puzzles: maybe instead of devising complex nonstandard ways of understanding the 
cosmological parameter values we should try to understand how the values of these 
cosmological parameters influence simple nonstandard dynamics. The relaxation 
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mechanism for the cosmological constant presented in this paper hopefully follows 
the latter route. 
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